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Abstract. This paper is devoted to the system of coupled KdV-like equations. It 
is shown that this apparently non-integrable system possesses an integrable reduction 
^-^ ■ which is closely related to the Volterra chain. This fact is used to construct the 

hyperelliptic solutions of the original system. 
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1. Introduction. 

The equation we are going to study is 

{dt - d^^^) Qn = Q {qn+iq'iqn-i) ^ ■ (1) 

The integrability/non-integrability of this apparently new system is now an open 
question. On one side, in this paper we estabhsh its relationship with the well-known 
integrable system, the Volterra hierarchy (VH) [H |2] (that is why we will call ([1]) 
the 'KdV- Volterra chain' (KdVVC)), and present the A^-phase periodic solutions from 
which one can derive the iV-soliton ones. According to the widely used hypothesis the 
existence of more than two-soliton solutions is a strong evidence of the integrability of 
the problem. On the other side, the system ([T]) can be easily reduced by taking g,„ = q 
to the generalized KdV equation, 

qt - qxxx = 30 g™g^. (2) 

with m = 4, that is not of solitonic type for m > 2 [3J. This statement does not have a 
rigorous proof and is based on the Painleve-like tests and numerical studies of interaction 
between solitary waves (see e.g. [1]). However, after so many years it seems improbable 
that the integrability of the generalized KdV equation can be established. Consequently, 
according to the Ablowitz-Ramani-Segur hypothesis (another integrability test that 
is not proved but has a long history of succesfull applications) the KdVVC is non- 
integrable. Here we have an apparent contradiction between the more than two-soliton 
and Ablowitz-Ramani-Segur tests: the former indicates the integrability of the KdVVC 
while the latter leads to the opposite conclusion. Probably we have a situation of 
the conditional integrability introduced by Dorizzi et al \E\ when equations in high 
dimensions (the higher KP equations in the case of [5]) are not individually integrable 
but become such, if one demands that their solutions also satisfy the lower equations of 
the hierarchy. As one can see below, our situation is almost the same. So, it is clear that 
the question of integrability /non-integrability of ([1]) is not trivial and deserves a special 
study, which is out of the scope of this paper. The main result of the presented work 
is the fact that equation ([1]) is one of a small number of nonlinear (2 + l)-dimensional 
systems for which an infinite family of explicit solutions can be derived. 

In this paper after reducing equation ([1]) in section [2] to VH (namely this reduction 
is the main topic of this paper) and discussing the hyperelliptic solutions for the latter 
(section [3]) we present the corresponding solutions for the KdVVC (sections H] and E]) . 

2. Reduction to the VH. 

The main result of this paper can be formulated as follows: if functions r„ solve the 

system 

[Dt — Dxxx) Tn ■ Tn-i = 3t„_|_2T„_3 — 6ar„+ir„_2 
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where a is a constant, D^^^ = D^ and Dt and D^ are Hirota's bilinear operators, 

df) [d^J n{t + t',x + x')v{t-t',x-x') ^ ^ ,(4) 
then the quantities 



Pn = -^ (5) 



r2 

' n 



solve the equation 

{dt -Cdcc- 9^xx) Pn = G {Pn+lplPn-l) ^ (6) 

with some function c = c{x,t). In the case of constant c (|6]) is related to ([1]) by simple 
Galilean transformation: g„(t, x) = Pn{t, x — ct). 

To prove this statement let us start from the first equation of system ([3]) which we 
rewrite as 

r^ Tn+l'T'n—2 /^\ 

OxfJ^n = (7) 

'n'n—1 

where 

/i„ = ln^^. (8) 

Tn-l 

From ([7]) we can obtain the following equations for the functions Pn dS]) that are related 
to fin by Pn = exp (/i„+i - /i„): 

^xPn = Pn (Un+l - M„) (9) 

dxxPn = Pn {Wn+1 - 4w„ + Wn-1 + M^+i + M^) (10) 

where, in order to make the formulae more readable, we have introduced 

(12) 
On the other hand, differentiating (jTj) we obtain 

dxx fJ-n = Wn - Wn-l (13) 

which implies 

dxxiTa.Tn = Wn + b (14) 

where b, in general, is a function of t and x but does not depend on the index n. 

Now we can calculate D^xx ^n ■ Tn-i and then dtfin and dtPn- Starting from the 
well-known formula for the Hirota's operators, 

^^^^=A... + A^ + 3A.(lna6),,, A = In ^, (15) 

after a little lengthy but simple calculations we get 

^xxx 'n ' 'n— 1 



n 




PnPn-1 




TnTn-l 










Jn 


= 


Un+lUn 


= 


Pn+lplPn- 


-1 


- In 


+2Tn- 
' n 


-2 



'n'n—1 



Un (Wn+l + 2Wn + 2Wn-l + Wn-2 (16) 

+M^+i + ul + u\_^ - 2m„,+iM„_i + 66) . 
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Substituting this result into the second equation of (^ we come to 

{dt - Cd^) fin = Un {Wn+l + 2Wn + 2Wn-l + Wn-2 (17) 

+M^+1 +ul + ul_^ + Un+lUn-l) 

with c = 6(6 — a), from which one can derive the expression for dt lnp„. However, before 
doing that it seems useful to present (ITTj) by means of the formulae for the derivatives 
of Un and Wn stemming from ([9]) 

dx Un = Un {Un+1 - Un-l) (18) 

dx Wn = Wn {Un+2 + Un+1 - Un - M„_i) (19) 

in the following two forms: 

(dt - C dx) /i„ = {dx +Un- Un-l) A„ + Bn-1 (20) 

= - {dx + Un+1 - Un) An + Bn (21) 

with 

An = Un {Un+1 + M„ + Un-l) (22) 

Bn = Un+lUn {2Un+2 + 3m„+i + 3u„ + 2m„_i) . (23) 

Using ( l20i) to calculate {dt — cdx) fin+i and ( l2Ti) to calculate {dt — cdx) fin we can easily 
obtain 

{dt - c 9^.) In pn = {dt-cdx) jin+i - {dt-c dx) fin (24) 

= {dx + Un+1 - Un) {An+1 + v4„) (25) 

or, recalling the fact that Un+i — Un~i = dx lnp„, 

{dt - C 9^.) p„ = dx Pn {An+l + An) . (26) 

Comparing the above expression with (TTOl) . 

P„ (^n+l + ^n) = <9xx Pn + 6w„p„, (27) 

we obtain 

{dt -cdx- dxxx) Pn = Qdx {WnPn) (28) 

which is nothing but (El). 

Equations ([3]) belong to the VH, one of the classical integrable models for which 
a wide range of solutions has already been constructed such as, e.g., solitons [U [2] and 
quasiperiodic solutions O [71 [H [9]. The main idea behind the present work is to use 
these results (with minor modifications) to derive solutions for ([1]) that is possible due 
to the established relation between ([3]) and ([T]). To illustrate this approach we obtain 
the periodic solutions for ([1]) starting from the hyperelliptic ones for the VH obtained in 
[7\ (see also [6] for the one-phase (elliptic) case). In the next section we rederive these 
solutions using an alternative to the methods in |71 El [9] and present the formulae we 
need for our purposes. 
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3. Hyperelliptic solutions for the VH. 

Consider tlie tiyperelliptic Riemann surface F 

29+1 

T: w' = V2,+i{0 = U ^^ - ^^^ (29) 

1=1 

wliicli is a compact Riemann surface of the genus g. One can choose a set of closed 

contours (cycles) {flj, &i}i=i,...,g with the intersection indices 

Qi o aj =hiO hj = 0, Oj o hj = 6ij i,j = l,...,g (30) 

and find g independent holomorphic differentials Uk satisfying the normalization 
conditions 



ujk = 5ik, i,k = l,...,g. (31) 

The matrix of the 6-periods, 

^ik = (p (^k (32) 



determines the so-called period lattice, Lq = {m + Qn, m, n G Z^}, the Jacobian of 
this surface Jac(r) = C^/Ln {2g torus) and the Abel mapping V — > Jac(r), 

P ^ I Lo (33) 

JPo 

where P is a point of F, P = (w, ,^), uj is the ^f-vector of the 1-forms, uj = (wi, . . . , cu^)^, 
and Pq is some fixed point of F. 

A central object of the theory of the compact Riemann surfaces is the 6'-function, 

6(0 = 0{C^), 

9{C)= Yl ^^P {^^ (^' ^^) + 2vrz (n, A } (34) 

TIGZ9 

where (n, ^) stands for ^f^i^iO? which is a quasiperiodic function on C^ 

e{C + n) =6(0 (35) 

9{C + ^n) = exp j-TTz (n, ftn) - 27ii (n, c)] (C) (36) 

for any n G Z^. 

The calculations presented below are based on the famous Fay's trisecant formula 
[Tot [TT] that can be written as 

4l 41 ^SS - 41 4: 0^1 e;^^ + e% 4: ^g ^S = 0. (3T) 

Here 

/ "^ rQ. \ 

(38) 





/ ."". Z"^' 


0^:3: 


= e c + l^/ ^ 




V ^=1 -^p^ 
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and the skew-symmetric function ep, ej, = —£q, is closely related to the prime form 
[TT] and is given by 



e^ = e{e+ u;\ (39) 

where e is a zero of the 6'-f unction: 6 {e) = 0. 

In what follows we will fix the points Pi, P2 and P3 as the images of the points 
^ = and ^ = 00, 

Pi = O = i+wo,0) 

P2 = 00 Wo = ^V2g+l{0) (40) 

P3 = = (-wo,0) 
and will consider the fourth point appearing in (13 7p as a variable one. Definition (HUj) 
leads to the possibility to take such integration paths in f l37|) that give 



00 00 



a;+ / a; = (41) 

o o 

(to this end it is enough, for example, to connect the points O and O with 00 by the 
curves in F having the same projections on the ^-plane). 

Now, by elementary transformations, we can convert the bilinear identity (|37|) into 
the form which gives us solutions for the VH. 

Making the shift C ~^ Cn = C + '^^ and introducing the function 0„(P) by 

0n(P)=^(C„ + '5(P)) (42) 

where the vectors u and S{P) are given by 
00 p 

iy= I uj, 6{P) = I iv (43) 

o o 

the Fay's identity ( 137|) can be rewritten as 

u{p) Qn-iio) e„(p) - e„(0) 0„-i(P) = u{P) e„_2(0) 0„+i(p) (44) 

with 

^(P) = ?!' § ^(P) = i- (45) 

This equation implies that the function 

UP) = a-'/'[u{P)rQr.{P) (46) 

where a is a constant that will be defined later solves 

T„_i(0) Tn{P) - Tn{0) T„_i(P) = UP) Tn-2{0) Tn+i{P) (47) 

with 



KdV-Volterra chain. 7 

Noting that the fraction Eq e^ / (e^) has the same zeroes and poles as the projection 

e(P): P = (^,0^e (49) 

one can conclude that it is possible to find the value of the constant a which ensures 
UP)=aP)- (50) 

Now we have to introduce the dependence of <^ on an infinite number of 'times' tj 
in such a way that the shift 

C^C + S{P) (51) 

is the Miwa's shift 



C(t) = C(...,t„...)^C(t + [e]) = C(...,t, + ^,...l. (52) 

To do this, let us take ^ as a local parameter near the point O of F. Thus the forms cu, 
can be presented as 

oo 

u, = u;,{Od^ = Y,''^ked^ (53) 

fc=0 

which leads to 

Ui= Ui{'q)dri = y^uJik—-— (54) 

Jo Jo j^^Q «+i 

and 

^(^) = ECA (55) 

where C,- is the vector with the components oJij-i- Now it is easy to check that if we 
take 

oo 

c(t) = E^.^^' (56) 

i=i 
then 

C(t + [e])-C(t) = 5(P). (57) 

In a similar way using the expansion 

lnu{P) = J2H' (58) 



(59) 







1 


one can 


introduce the function 




0(t) 


j 


with 




J 




0,= 


J\ 



(60) 
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such that 

u{P) = e'^(t+Kl)-*W. (61) 

Combining the above results one can present the functions Tn{P) as functions r„ (t): 

T„(0) = r„(t), T„(P) = r„(t + [^]) (62) 

where 

r„ (t) = a'^'/2gn0(t)^ f. ^^) ^ ^^j ^g3) 

and conclude that functions r„ (t) solve 

Tn-i (t) r„, (t + [e]) - r„ (t) r„_i (t + [^]) = ^ r„_2 (t) r„+i (t + [^D . (64) 

Equation flM|) is nothing but the so-called functional representation of the VH [12j. 
Expanding it in the power series in ^ one can show that any of its solutions also solves 
equations of the VH, the first three of which are given by 

Ai Tn ■ r„_i = r„+ir„_2 (65) 

Dt2 Tn ■ Tn-l = Dti T^+i ■ r„_2 (66) 

(8 A3 + Dt^hh) Tn ■ Tn^l = {GDt^ + SDt^ti) ^«+l ' Tn-2 (67) 

Eliminating the t2-derivatives one can rewrite the last equation as 

(A3 - Aititi) Tn ■ Tn-1 = 3Tn+2Tn-3 " 66r„+ir„_2 (68) 

where, recall, b = dt^ti lnT„ — Wn and does not depend on n by virtue of (165|) . Thus the 
functions r„ given by (!63l) together with (J56l) and (J58ll . ( !59l) are hyperelliptic solutions 
for equations ([3]) with a = b. This means that we have everything necessary to construct 
the hyperelhptic solutions for the KdVVC ([1]). 

4. Hyperelliptic solutions for the KdVVC. 

As it can be seen from the results of the previous section, the r-functions ( l63ll 
constructed from the hyperelliptic ^-functions satisfy equations ([3]), (fl^ with a = b. 
This means that we do not need to know the exact values of the constants a and b 
because the parameter c appearing in IQ for pn is automatically equal to zero in this 
case. In other words, we do not need to make the Galilean transformation and the 
hyperelliptic solutions for the KdVVC, g„, are given by 

qn = Il^±l^. (69) 

T 
n 

Recalling the definition of r„ 063 p one can note that the functions (t) cancel 
themselves due to the structure of fl69|) . We have introduced the function u{P) in 
fH6|l . and hence the function (p (t) in fl6T|) . to rewrite the Fay's identity as the standard 
bilinear functional equations of the VH p!2]. In principle one can deduce directly 
from (Hi|) some functional equations for the quantities 9„_i9„+i/0^ which lead to the 
differential equations solved by Pn and g„. However, in this case one has to perform more 
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cumbersome calculations giving up the advatages of the bilinear approach. Moreover, 
we have not used this way because one of the goals of this work is to demostrate the 
relationship of ([T]) with the standard Volterra equations. The disappearance of the factor 
£"</>(*) from the final formulae for g„ is a typical effect of the bilinear and the inverse 
scattering approaches: periodic solutions of integrable equations lead to quasiperiodic 
solutions of corresponding linear and bilinear problems. In some sense, the Fay-like 
formulae of the previous section contain more information than we need for the present 
work: for example, the above mentioned function (t) which is inessential for g„ would 
be crucial if we were deriving the Baker- Akhiezer function related to (!69|) . 

The role of the parameter a is different: the factor a" /^ 'survives' the bili- 
nearization procedure of calculating Tn-iTn+i/r^ and appears in the final formulae as 
the 'amplitude' of the solutions. 

Thus, 

^ (C+l) (Cn-l) . . , .^^x 

(C„) 
The dependence of g„ on x and t can be established by comparing (^ and (^E\i . It is 
obvious that we have to identify x with ti and t with ^3 neglecting all the other times 
tk, k ^ 1,3 (they play the role of a constant): 

C{x,t) = CiX + Cs^ + C*, C* = constant (71) 

where (^^ 3 are the vectors with the components 

(Ci)i = ^io = ^i{0) 

(C3). = ^^2 = Wm- 

5. Elliptic solutions for the KdVVC. 

In this section we rewrite solutions ( ITOl) . ( 17T]) and calculate all the constant parameters 
in terms of the Jacobi elliptic functions. Expression ( !70|) for g = 1 becomes 

^3 (Cn+l) ^3 (Cn-l) ^ ^ , .„„s 

qn = a -J-— , („ = ( + nz/ (73) 

V3 yi^n) 

(we use the notation of the book [13]) and can be rewritten using the Jacobi functions 
as 

qn = q*fizn) (74) 

where 

f{z) = f{z; a,k) = l-k^ sn\a, k) sn^{z, k) (75) 

and 

Zn = K{k) (2C„ + 1) , a = 2K{k)u (76) 
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with K{k) being the complete eUiptic integral of the first kind. Starting from the 
standard formulae for the derivatives of the elliptic functions one can obtain the following 
identity for the function f{z): 

f"iz) = P,fiz)+PJ{z) + (3o (77) 

where 

B - 6 

sn a 

Pi = ^-4(1 + A;2) (78) 

sn a 

(3o = ^ + A(l + k^)-2esn^a. 

sn a 

Another useful property of the function f{z) stems from the superposition formulae for 

the elliptic functions: 

f{z + a)f{z)f{z -a)= ^,f{z) + 70 (79) 

where 

71 = 4 cn^ a dn^ a 

(80) 
7o = P{a) — Acn^ adn^ a. 

Using (1771) and (179|) one can verify the fact that the functions 

qn{x,t) = q^f {z + na) (81) 

with 

z = z{x,t) = Kx + zut + constant. (82) 

satisfy the equation 

K^—r-T - ^qn + 6fi;g„_ig^g„+i = C = constant (83) 

provided 



k' = -6^gt w = Pii^' + Q-fo^qt (84) 

P2 

(the constant C is given by C = 7oK^g*) which implies that functions g^ solve ((1]). 
Substituting the values of /3j, 7^ we come to the final result: the elliptic solutions of the 
KdVVC are given by dH]), ((75]) and (182]) with 

K, = 2 sn a cn a dn a g^ (85) 



-cu = 6 



1 12, ,., 



sn^ 2a sn^ a 3 



K^ (86) 
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